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1 Introduction 

Infrared radiance measurements were acquired from a radiometer on the NASA ER-2 during a 
coincident Landsat 5 overpass on October 28, 1986 as part of the FIRE Cirrus IFO in the vicinity of 
Lake Michigan. A comparative study is made to infer microphysical properties of the cirrus cloud 
field. Radiances are derived from the image by convolving the ER-2 radiometer’s effective field 
of view along the flight path. A multistream radiative transfer model is used to account for the 
differences in spectral bandwidths, 10.40-12.50 /xm for the Landsat band and 9.90-10.87 /xm for the 
radiometer. 

2 Instruments 

The primary aircraft based instrument employed for this study is a two spectral channel nar- 
row field of view radiometer (NFOV). This instrument was mounted within a pod attached to the 
fuselage of the NASA ER-2 aircraft. Both channels detect upwelling infrared radiation within a 
conical field half-angle of 8.3 degrees. The optically sensitive components are electrically calibrated 
pyroelectric detectors (Valero et al. 1982; Geist and Blevin 1973). A reflective chopper operating 
at 18 hz alternately exposes the detectors to external radiation and a liquid nitrogen cooled zero- 
radiation reference. The temperature corrected spectral bandwidths (full width at half maximum) 
of the interference filters are 6.14-7.14 and 9.90-10.87 /xm. In this study only the 10 /x in channel 
bandwidth overlapped one of the Landsat 5 bands: TM-6. 

The basic operation principles of the Landsat thematic mapper are described elsewhere (Engel 
and Weinstein 1983). Standard relations were used for converting the digital pixel counts to radi- 
ance and brightness temperature (Markham and Barker 1986). The Landsat 5 TM6 band has an 
equivalent rectangular spectral bandpass of 10.40-12.50 /x m (Wielicki and Parker 1987). Radiances 
derived from the Planck function integrated over this bandpass agree very well with the Markham 
and Barker formulas for which this integral is approximated by a simple analytic expression. 

3 Analysis 

In order to make meaningful comparisons between our NFOV ER-2 radiometer and Landsat 
image data, a systematic procedure was used to prepare both data sets. First, a geometrically 
registered satellite image was obtained. Next, the aircraft flight tracks derived from navigational 
data were adjusted for cloud level winds. Radiances were determined from the image by convolving 
the ER-2 radiometer’s effective field of view with image pixel values along the flight track. The 
previous two steps require an estimate of the altitude of the optically dominant cirrus layer; this was 
aided by lidar data. The resulting signals were expressed in spectral bandpass independent units of 
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brightness temperature verses aircraft coordinates. Finally the two sets of radiance measurements 
were prepared for a 2 channel radiance correlation plot for comparison to model calculations. 

4 Model 

A discrete ordinates radiative transfer code designed for multi-layered plane-parallel media 
(Stamne8 et al., 1988) was used to calculate radiance values at the corresponding altitudes of the 
measurements, 19.1 km for the ER-2 and >70 km for Landsat. The number of computational polar 
angles (streams) was set to 8, as larger values did not generally yield significant differences in the 
results. 17 atmospheric layers were chosen so as to achieve reasonable homogeneity in temperature 
(less than 10°K variation) and particle composition within each layer. The upwelling radiance was 
evaluated at a polar angle of 0 degrees to correspond to the conditions of the measurements. 

The assumed ice particle size distributions are described in terms of Mie spheres of a single 
equivalent radius. Mie calculations were done for a selection of water and ice sphere radii using ex- 
perimentally determined values for the complex index of refraction of ice (Warren, 1984) and water 
(Downing and Williams, 1975). Particle concentrations were chosen to yield the desired zero zenith 
angle optical depths for extinction through all the cloud layers at a reference wavelength of 11.4 pm. 

The atmospheric profile input to the model was derived from the 1500 GMT Greenbay radiosonde. 
Since direct measurements of the Lake Michigan surface temperature were not available the Green- 
bay surface temperature was used: T, = 281.9°K. The lake surface was assumed to be Lambertian 
with a wavelength independent albedo of 0.04. ER-2 based lidar measurements (Spinhirne et al. 
1988) indicate that there are two distinct cloud layers over most of the Landsat image extract, with 
altitude ranges of about 7 to 8 km and 9 to 11 km. 

5 Results and Conclusions 

Fig. la is the Landsat image extract. Fig. lb shows the ER-2 radiometer and Landsat signals 
along the flight track in units of brightness temperature. Figs. 2 and 3 are correlation plots. The 
ER-2 and satellite scales are expressed in radiance units (Wm^sr -1 ) for the instruments’ respective 
bandpasses. The measured data points of Fig. 2 (represented as small dots in Fig. 3) are derived 

from the data represented in Fig. lb. Data is distinguished as to ice or water according to the lidar 
plot of Fig. lc. 

The model calculations are plotted as connected sets of specific symbol types. Each set cor- 
responds to radiance results for a range of optical depths using particles with particular radius 
selections. Beginning at the upper right convergence point the optical depth values (referenced at 
11.4 pm) are 0, 0.125, 0.25, 0.375, 0.5, 0.625, 0.75, 1.0, 1.25, 1.5, and 2.0. For the standard model 
(connected unfilled symbols in Fig. 2 and solid symbols in Fig. 3) a two layer cloud profile was 
chosen with ice spheres in the upper layer having about half the effective radius of spheres in the 
lower layer. The 11.4 pm optical depths of the lower and upper layers were set equal. Unfilled 
symbols in Fig. 3 represent a profile with water drops present in the lower layer and nothing in 
the upper. The following conclusions are based on comparisons of measured radiances with model 
results on the correlation plots: 

1. Assuming the above standard model, ice particles with equivalent sphere radii of 8 to 15 fx m 
for the upper layer and 15 to 30 pm for the lower layer may be inferred. This contrasts with particle 
sizes typically greater than 100 pm determined by in-situ measurements. The lack of sensitivity 
of the in-situ probes to small ice crystals and the uncertainties in ice-water content determinations 
causes difficulties in resolving whether a large number of small particles are actually present as 
opposed to complex large particles with small protrusions. The combined cloud extinction optical 
depth for both layers at the reference wavelength of 11.4 pm ranges from about 0.3 to 2. 

2. Assuming a 1-layer lower water cloud yields results consistent with approximately 8 pm radii 
water droplets. This compares with in-situ measurements of radii about 4 pm. Distinguishing mea- 
surements made with water clouds present from those made for just ice clouds gives no apparent 
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separation on the correlation plots. Thus water and ice clouds cannot be distinguished from the 
correlation plots alone under the conditions of this case study. The optical depth scale is only mod- 
erately changed (< 20%) from the 2-layer ice model results. 

3. Limits on the estimated instrument error may be established by means of self-consistency of 
the model results. Specifically, if Landsat radiances values were too high, the inferred water cloud 
droplet sizes would be unrealistically large. Thus our retrieval of relatively small ice crystals cannot 
be attributed to a Landsat (nor an ER-2 radiometer) calibration error. 
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